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Abstract. A combinatorial substitution is a map over tilings which allows to 
define sets of tilings with a strong hierarchical structure. In this paper, we show 
that such sets of tilings are sofic, that is, can be enforced by finitely many local 
constraints. This extends some similar previous results (Mozes'90, Goodman- 
Strauss'98) in a much shorter presentation. 



1. Introduction 

Tiling some space with geometrical shapes, or tiles, consists into covering this 
space with copies of the tiles. When a set of tilings can be characterized by adding 
finitely many local constraints on tiles so that the set of tilings corresponds exactly 
to the set of tilings satisfying the constraints, such a set of tilings is called sofic. 
Soficity corresponds to the interesting idea that the validity of a tiling can be locally 
proved by decorating tiles with constraints. In this paper, we contribute to a general 
question: which sets of tilings are sofic? 

Substitutions provide a simple way to express how a set of tilings can be ob- 
tained by iteratively constructing bigger and bigger aggregates of tiles. The strong 
hierarchical structure of substitutions tilings permits to enforce global properties, 
for example aperiodicity. To prove that the set of tilings generated by a substitution 
is sofic, one has to encode the global hierarchical structure into local constraints on 
tiles. Such technique is at the root of classical papers on the undecidability of the 
Domino Problem [21 E] . In the case of tilings on the square grid, Mozes [6] proved 
in a seminal paper that the set of tilings generated by rectangular non-deterministic 
substitutions satisfying a particular property is sofic. Goodman-Strauss [5] proved 
that such a construction method can be extended to a wide variety of geometrical 
substitutions. In this paper, based on ideas developed in [7], we further extend the 
construction to a broader class of substitutions by replacing the geometrical condi- 
tions by combinatorial conditions while decreasing the length of the presentation. 

Let us sketch some definitions to state our main theorem. A sofic tiling is a 
valid tiling by a finite set of tiles with decorations. Combinatorial substitutions, 
introduced by Priebe- Frank in [8], map a tiling by tiles onto a tiling by so-called 
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macro-tiles (finite tilings, here assumed to be connected), so that the macro-tiles of 
the latter are arranged as the tiles of the former. The limit set of a substitution is 
the set of complete tilings that admits preimages of any depth by the substitution. A 
good substitution is a substitution that is both connecting (a combinatorial condition 
ensuring that there is enough room for all information to flow) and consistent (a 
geometrical condition ensuring that the substitution can be correctly iterated). The 
main result obtained in this paper is: 

Theorem 1.1. The limit set of a good combinatorial substitution is sofic. 

The paper is organized as follows. Sections [2] and [3] formally define main no- 
tions, in particular sofic tilings and good combinatorial substitutions. Section [4] 
then presents self-simulation, which plays a central role in the constructive proof of 
Theorem 1.1 which is given in section [5] Last, section [6] concludes the paper by 



discussing an important parameter of this proof. 

A single example illustrates definitions and results throughout the whole paper: 
the "Rauzy" example. It relies on the theory of generalized substitutions introduced 
in pQ, a self-contained presentation of which is beyond the scope of this paper. Let 
us just mention that Rauzy tilings are digitizations of the planes of the Euclidean 
space with a specific given irrational normal vector. More details, as well as the 
results we here implicitly rely on, can be found in [31 H] - We chose this example, not 
the simplest one, because it is not covered by results in [6l [5]. 



2. Sofic tilings 

Polytopes are assumed to be homeomorphic to closed balls of lR d and to have 
finitely many faces, with the (d — l)-dimensional faces being called facets. 

A tile T is a polytope of the Euclidean space M. d . A tiling Q of a domain DcK d 
is a covering of D by interior-disjoint tiles, with the additional condition that two 
tiles can intersect (if they do) only along entire faces. A tiling is said to be finite if 
its domain is bounded. 

Two tiles are said to be adjacent if they intersect along at least one facet, and a 
tiling is said to be connected if any two of its tiles can be connected by a sequence 
of adjacent tiles. 

A facet of a tiling is said to be external if it is on the boundary of the domain, 
internal otherwise. We denote by dQ the set of external facets of a tiling Q. If this 
set is empty, then the tiling is said to be complete: its domain is the whole M d . 

A decorated tile T is a tile with a real map defined on its boundaries, called 
decoration. Two adjacent decorated tiles are said to match if their decorations are 
equal in any point of their intersecting facets. A decorated tiling Q is a tiling by 
decorated tiles which pairwise match (when adjacent). Decorations are thus local 
constraints on the way tiles can be arranged in a tiling. 

A tileset r is a set of decorated tiles. It is said to be finite if it contains only a 
finite number of tiles up to direct isometries. A decorated tiling whose tiles belong 
to a tileset r is called a r-tiling; the set of r-tilings is denoted by A T . 

One says that a tiling can be seen as a decorated tiling if both are equal up to 
decorations. One denotes by % the map which removes the decorations. One easily 
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checks that any tiling can be seen as a r-tiling if r can be infinite. The interesting 
case is the one of finite tilesets: 

Definition 2.1. A set of tilings is said to be sofic if it can be seen as the set of 
r-tilings of some finite tileset r. 

For example, one can wonder whether the Rauzy tilings are sofic (Fig. [I]). 



Figure 1: A Rauzy tiling (partial view). Are Rauzy tilings sofic? 



3. Combinatorial substitutions 

Definition 3.1. A combinatorial substitution is a finite set of rules (P, Q,^), where 
P is a tile, Q is a finite connected tiling, and 7 : dP — > dQ maps distinct facets on 
disjoint sets of facets. The tiling Q is called a macro-tile, and if / is the k-th facet 
of P, then 7(/) is called the k-th macro-facet of Q. 

Fig. [2] illustrates this definition. 

We call tiling by macro-tiles a tiling whose tiles can be partitioned into macro- 
tiles, with each macro-facet belonging to the intersection of exactly two macro-tiles. 
We can associate with each combinatorial substitution a binary relation over tilings: 

Definition 3.2. Let a be a combinatorial substitution. A tiling T by tiles of a 
and a tiling T' by macro-tiles of a are said to be a-related if there is a one-to-one 
correspondence between the tiles of T and the macro-tiles of T' which preserves 
the combinatorial structure, that is, such that the a-th facet of a first tile of T 
matches the 6-th facet of a second tile of T if and only if the a-th macro-facet of 
the first corresponding macro-tile of T' matches the 6-th macro-facet of the second 
corresponding macro-tile of T' . One calls T a preimage of T' and T' an image of T. 

For example, a Rauzy tiling can be uniquely seen as a tiling by Rauzy macro- 
tiles (Fig. [3]) and has a unique preimage, which turns out to be itself a Rauzy tiling 
(both facts are non-trivial; they follow from results proven in [3]). 

In particular, the relations associated with combinatorial substitutions yield a 
strong hierarchical structure on so-called limit-sets: 
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Figure 2: These five rules define the so-called Rauzy combinatorial substitution (a 
facet / and the corresponding macro-facet 7(/) are similarly marked). 




Figure 3: A Rauzy tiling can be uniquely seen as a tiling by Rauzy macro-tiles. 

Definition 3.3. The limit set of a combinatorial substitution a, denoted by A CT , is 
the set of complete tilings which admit an infinite sequence of preimages. 

For example, the limit set of the Rauzy combinatorial substitution is exactly the 
set of Rauzy tilings (again, this non-trivial fact follows from results proven in [3]). 



Let us now turn to the good combinatorial substitutions to which Theorem 1.1 
applies. First, a good combinatorial substitution must be connecting: 

Definition 3.4. A combinatorial substitution o is connecting if, for each rule 
(P,Q,j), the dual grapb[]of Q has a subgraph N, called its network, such that 
(1) A is a star with one branch for each macro-facet, and the leaf of the k-th 
branch is a tile with a facet, called k-th port, in the k-th macro-facet of Q; 

lr The dual graph of a tiling is the graph whose vertices correspond to tiles of the tiling and 
whose edges connect vertices corresponding to adjacent tiles. 
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(2) Each macro-facet has non-port facets, and removing the edges of N and its 
central vertex yields a connected graph which connects^] all these facets; 

(3) the center of N corresponds to a tile in the interior of Q, called central tile; 

(4) whenever two macro-tiles match along a port, they also match along the 
corresponding macro-facet. 

Informally, the two first conditions ensure that the macro-facets are big enough 
to transfer via the network all the informations (encoded by decorations) that we 
need to enforce the hierarchical structure of the limit set. In particular, one easily 
sees that connectivity could not be achieved for combinatorial substitutions on the 



real line: Theorem IT can apply only in dimension two or more. The third condition 
ensures that, by iteratively considering macro-tiles of macro-tiles (that is, when going 
higher and higher in the hierarchy of the limit set), we get tilings covering arbitrarily 
big balls. The last condition associates a port with its macro-facet (it is equivalent 
to the "sibling-edge-to-edge" condition of [3]). 

Second, a good combinatorial substitution must be consistent: 

Definition 3.5. A combinatorial substitution a is said to be consistent if any tiling 
by macro-tiles of a admits a preimage under a. 

Intuitively, consistency ensures that, if a tiling meets all the combinatorial con- 
ditions to have a preimage, then there is no geometrical obstruction to the existence 
of such a preimage. 

Open Problem 3.6. Characterize consistent combinatorial substitutions. 

For example, the Rauzy combinatorial substitution is connecting (one easily 
finds a suitable network for each rule, see Fig. [4] and [5]) and consistent (this non- 
trivial fact follows from results in |4], where explicit maps defined over tilings are 
shown to be equivalent to such combinatorial substitutions). Theorem thus 
yields that its limit set, i.e., the set of Rauzy tilings, is sofic. 






Figure 4: The Rauzy combinatorial substitution is connecting. 



One says that a subgraph connects a set of facets if these facets all belong to tiles which 
correspond to vertices of this subgraph. 
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Figure 5: A Rauzy macro-tile (top-left) and an image of it under the Rauzy combina- 
torial substitution, that one could call Rauzy "macro- macro-tile" (center). 

4. Self-simulation 

Definition 4.1. Let a be a combinatorial substitution with the rules {(Pi, Qi, 7i)}j- 
A tileset r is said to a -self -simulates if there is a set of r-tilings, called r-macro-tiles, 
and a map <p from these r-macro-tiles into r such that 
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(1) for any r-macro-tile Q, there is i such that tt(Q) = Qi and vr(0(Q)) = P^ 

(2) any complete r-tiling is also a tiling by r-macro-tiles; 

(3) the a-th macro-facet of a r-macro-tile Q can match the 6-th macro-facet of 
a r-macro-tile Q' if and only if the a-th facet of the r-tile <fi(Q) can match 
the 6-th facet of the r-tile 4>(Q')- 

Proposition 4.2. If a tileset a -self- simulates for a consistent combinatorial substi- 
tution a, then its complete tilings are, up to decorations, in the limit set of a. 

Proof. Consider a complete r-tiling V. Conditions (l)-(2) ensure that removing 
the decorations of V yields a tiling by macro-tiles of a, say P. The consistency 
of cr ensures that P admits a preimage under a, say R. Let us show that R can 
be endowed by decorations to get a r-tiling. Consider a tile T of R. This tile 



corresponds (via the one-to-one correspondence in Def. 3.2) to a macro-tile Q of P, 
which is itself the image under tc of a r-macro-tile Q of V. We associate with T the 
r-tile 0(<2). Condition (3) ensures that replacing tiles in R by their such associated 
r-tiles yields a complete r-tiling, say 1Z. We can repeat all this process, with 1Z 
instead of V . By induction, we get an infinite sequence of complete tilings, with 
each one being the preimage under a of the previous one. Thus, V G A CT . ■ 



5. Constructive proof of Theorem 1.1 



Let a be a good combinatorial substitution with rules (Pi, Qi, 7^. We here rely 
on the fact that a is connecting to construct a finite tileset r which a-self-simulates. 
The consistency of o then ensures, via Prop. 4.2, that 7r(A r ) C A CT holds. We also 



show that the converse inclusion holds. This thus constructively proves Theorem 1.1 



5.1. Settings 

Let Ti, . . . , T n and fi, ■ ■ ■ , f m be numberings of, respectively, the tiles and the 
internal facets of all the Qj's. Given the k-th facet of the tile T i; N a [i,k) stands 
either for its index if it is an internal facet, or for a special value "port", "macro- 
facet" or "boundary" otherwise (depending whether it is a port, a non-port facet in 
a macro-facet or another external facet). 

Each tile of r is a Tj endowed with a decoration which encodes on each facet 
a triple (f,j,g), where / and g are either facet indices or special values "port", 
"macro-facet" or "boundary", and j is either zero or a tile index. We call / the 
macro-index, j the parent-index and g the neighbor-index. This clearly allows only 
a finite number of different tiles. 

We skip the technical details concerning the way these triples are encoded by 
decorations^} We assume that two decorated tiles match along a facet if and only if 
the same triple is encoded on both facets, and that the only direct isometry which 
leaves invariant the decoration of a facet is the identity, so that a decorated tile 
cannot trivially match with a translated or rotated copy of itself. 



'Recall that the decoration of a tile is a real map defined on its boundary. 
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5.2. Decorations 

The five following steps completely define a tileset r. 

1. The macro-index of the k-th facet of any decorated Tj is N a (j,k). This step 
ensures that any complete r-tiling can be uniquely seen as a tiling by r-macro-tiles. 

2. Consider a decorated non-central tile of Qi. Its facets which are internal and not 
crossed by the network have all the same parent-index, also called parent-index of 
the tile, which can be any j such that Tj = Pj. Its facets which are external, port 
excluded, have parent-index 0. This step ensures that the r-tiles of a r-macro-tile 
Q share a common parent-index j; the tile Tj is called the parent-tile of Q. 

3. Consider, in a non-central r-tile with parent-index j, a facet which is neither a 
port nor crossed by the network. Its neighbor-index is either N a (j, k) if it is in the 
k-th macro-facet, or equal to its macro-index otherwise. This step ensures that the 
macro-facets (ports excepted) of a r-macro-tile are equivalent to the macro-indices 
of its parent tile (once decorated). 

4. Consider a non-central r-tile with parent-index j. Its facets which are either 
A;-th port or crossed by the k-th branch of the network have all the same pair of 
parent/neighbor indices: it can be any pair not forbidden on the A;-th facet of Tj 
(Steps 1-3). This step ensures that each port of a r-macro-tile is equivalent to the 
pair of parent /neighbor indices of a decorated parent-tik^J 

5. Whenever a non-central r-tile T has as many facets as a central tile Tj, we define 
a central r-tile T' by endowing each k-th facet of Tj with the pair of parent/neighbor 
indices on the k-th facet of T (the macro-indices are defined as usual, see Step 1). 
One says that T' derives from T. 



5.3. First inclusion 

Let us show that the above defined tileset a-self-simulates. Given a r-macro-tile 
Q with parent-index j and central r-tile T 7 , let 0(<2) be the decorated tile obtained 
by endowing the k-th facet of Tj with the parent /neighbor indices on the k-th facet 
of T' (the macro- indices are defined as for any tile, see Step 1). One checks that 4> is 



a map satisfying conditions (l)-(3) of Def. 4.1 It remains to check that <fi(Q) G r. 

Let T be a non-central r-tile from which derives T' ■ If Tj is central, then <p(Q) 
also derives from T, hence is in r. Otherwise, Step 4 ensures, for each k, that 
the parent /neighbor indices on the k-th facet of T' appear on the k-th facet of a 
decorated Tj. In particular, this holds on facets of Tj which are not crossed by a 
network. Since the neighbor-indices of these facets can onlj|^] appear on a Tj (see 
Steps 1 and 3), T is a decorated Tj. This yields 0(<2) = T, and thus 0(Q) is in r. 



Prop. 4.2 then applies and yields the first inclusion 7r(A r ) C A c 



With the problem that a parent-tile can be decorated in different ways, and nothing yet prevents 
the ports from mixing these decorations. 

5 Actually, a facet-index appears on the two tiles of a macro-tile which share the corresponding 
facet. We thus need, in order to completely characterize Tj, either to assume that there is at least 
two such facets (this is a rather mild assumption), or to endow facets with an orientation and to 
allow two facets to match if and only if they have opposite orientations. 
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5.4. Second inclusion 

Let us extend r in a tileset r' as follow. For each r-tile T and each subset S of 
its facets, we define a r'-tile Ts by replacing the decorations of the facets in S by a 
special decoration "undefined". 

Now, let P be a tiling in A CT . Consider an infinite sequence (P n ) n >o of successive 
preimages of P. Given n > 0, P n can be seen as a r'-tiling: it suffices to endow any 
facet with "undefined". Then, P n -\ can be seen as a tiling by r'-macro-tiles, with 
"undefined" decorations appearing only on the network. Indeed, consider a macro- 



tile of P n -\ which corresponds (via the one-to-one correspondence in Def. 3.2) to 
a tile Tj in P n : it suffices to endow its tiles as in steps 1-3 of the definition of r, 
with Tj being the parent-tile, and with the decoration "undefined" on the facets 
crossed by the network. This can be iterated up to P = Pq, and the third condition 



of Def. |3.4| ensures that the decorations "undefined" appear only on sort of grids 
whose cells have bigger and bigger size. 

Thus, by making n tend to infinity, one can see P as a r'-tiling whose "undefined" 
decorations, if any, form either a star with k infinite branches, or a single biinfinite 
branch. In the first case, we can replace the central r'-tile by any r-tile with A; 
facets, and then the tiles on branches by r-tiles which carry decorations to infinity. 
In the second case, we can replace the r'-tiles by r-tiles which carry any decoration 
on the whole branch. In any case, we can thus see P tiling. 

We thus have the second inclusion A CT C 7r(A T ). Both inclusions prove Theorem 

6. On the number of tiles 

Let us conclude this paper by discussing the size #r of the tileset r defined in 



the previous section (although its finiteness suffices for Theorem 1.1). 



Consider a good combinatorial substitution with r rules. Fix a network as in 



Def. 3.4 for each of these rules. Let n and m denote the total number of, respectively, 
tiles and internal facets of the macro-tiles of these r rules. Among these n tiles, let 
p denote the number of those which lie on a network. 

First, the n — p tiles not on the network can be decorated in at most n different 
ways, according to the parent-index they carry. This yields at most Nq = (n — p)n 
r-tiles. Then, each of the p — r non-central tiles on the network can be decorated in 
at most (n—p)n+pnm different ways: n — p parent-indices which correspond to tiles 
not on the network, hence allow at most n pairs of parent /neighbor indices carried on 
the network, and p parent-indices which correspond to tiles on the network, hence 
allow at most nm pairs of parent/neighbor indices carried on the network. This 
yields at most N p = (p — r)((n — p)n +pmn) r-tiles. Last, the r central tiles can be 
decorated in at most as many different ways as there is non-central r-tiles. Finally, 
this yields 

#r < (r + 1)(JV + N p ) < (r + 2)p 2 mn. 
This bound is, for example, about one billion in the Rauzy case. 



In order to reduce this huge number of tiles, it is worth noting that, instead 
of carrying a parent-index through all the tiles of a macro-tile, it suffices to carry 
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Figure 6: A second network for the Rauzy combinatorial substitution. 

it along a second network connecting the macro-facets of this macro-tile and inter- 
secting each of the branches of its (first) network (see, e.g., Fig. [6]). The "control" 
described in Step 4 is then performed only on tiles where both networks crosses, 
while we simply allow any possible pairs of parent/neighbor indices to be carried 
through the tiles which are only on the first network. If we denote by q the to- 
tal number of tiles on these new second networks and by c be the total number of 
crossings between second and first networks, a similar analysis yields 

#r < (r + 1)(A^ + N' q + N' p + N' c ) < (r + 2)cp(m + qn), 

where: 

N' Q = n - p - q + c, N' p = (p - c)(m + qn), 

K = (q-c)n, N' c = c(N^ + N' q + N' p ). 

This bound is, for example, about 70 millions tiles in the Rauzy case. 

This last bound is huge but generic. One can hope to dramatically decrease 
this bound in specific cases. Indeed, most of the r-tiles correspond to tiles which 
simply carry any possible information (the tiles on networks, crossings excepted). 
Since these tiles all play the same role, it would be worth to replace all the tiles 
on a network by a single tile (one can thus hope to gain a factor pq in the above 
bound - this would yield about 25000 tiles in the Rauzy case). This shall however 
be done carefully, so that tiles still necessarily form macro-tiles. Note that it should 
be much easier in dimension d > 3, since the cohesion of macro-tiles can be more 
easily enforced without relying on tiles on networks. 
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Appendix A. Complete example: the 3x3 square substitution 
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Figure 7: a simple substitution 
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Figure 9: chosen network and numberings 
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Figure 10: decoration of facets 
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Figure 11: initial decoration scheme 
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Figure 12: after step 1: macro-indices are fixed 
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ure 13: after step 2: parent-indices are fixed outside network 
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where n = N a {j, N), s = N a {j, S),w = N a (j, W) and e = N a {j, E). 
Figure 14: after step 3: all decorations fixed outside network 
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parent /neighbor indices 
Figure 15: after step 4: all decorations fixed but for central tiles 
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where, for any j, n = N a (j,N), s = N a (j,S), w = N a {j,W), 
e = N a (j,E), (cej,f3j) are valid parent/neighbor indices for a tile 
Tfc with parent-index j 



Figure 16: after step 5: all decorations are fixed 



